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Abstract 

We studied the dynamics of surfacea and wake charges induced by a charged particle travers- 
ing a dielectric slab. It is shown that after the crossing of the slab first boundary, the induced 
on the slab surface charge (image charge) is transformed into the wake charge, which overflows 
to the second boundary when the particle crosses it. It is also shown, that the polarization of 
the slab is of an oscillatory nature, and the net induced charge in a slab remains zero at all 
stages of the motion. 
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1 Introduction 



As it passes through a medium, a fast charged particle excites oscillations of the charge density 
behind itself [1-3] . These wakefields and the particle energy losses associated with their excita- 
tion have been studied widely for a variety of media [3-8]. Wakefields have recently reattracted 
interest because of the development of new methods for accelerating particles [9-11]. 

In most studies of wakefields it has been assumed that the medium is unbounded. The wake- 
fields are excited as the particle enters the medium, or they disappear when the particle leaves 
the medium, because of various transient polarization processes which occur near the inter- 
face. Among these processes, the excitation of surface oscillations and the associated additional 
energy loss have been studied previously [12-20]. In connection with the development of new 
particle acceleration methods, numerical calculations have determined the distance from the 
sharp plasma boundary at which the amplitude of the wakefield excited by an ultrarelativistic 
particle reaches the same level as in an unbounded medium [21]. 

Fairly recently, in connection with problems of emission electronics and optoelectronics, the 
image charge and the dynamical image potential created by a moving particle has also been 
investigated. To describe the process of formation of the image charge, various approaches 
(quantum mechanical, the hydrodynamic, etc.) and various models of the medium have been 
employed [22-25]. 

In the present paper we analyze the dynamics of reversal of the sign of the charges induced 
at the slab boundary (repolarization of the slab) as the particle crosses the interface. The 
process is found to be of a nonmonotonic, oscillatory nature. The case of normal incidence of 
a particle through the slab is considered. 

The paper outline is as follows. In section 2, general expressions for the density of wake 
charge and total wake charge have been found, using Poisson's equation. In section 3, general 
expressions for the density of induced surface charges and total charges have been found, using 
expressions for the normal component of the electric field in the internal and external space of 
the slab [26] . We apply the results obtained in sections 2 and 3 to the case when slab constructed 
from a diatomic cubic ionic crystal or polar semiconductor. In section 4, the obtained results 
are discussed. 
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2 The electromagnetic field of charged particle traversing a slab 

We consider a fast particle of charge q moving with a velocity u along the z-axis normal to the 
boundaries of a slab characterized by a local dielectric function e(u). The time interval t during 
which the particle moves through the medium is < t < a/u, where a is the slab thickness. 
Outside this interval the particle moves in a vacuum. 

Ginzburg and Tsytovich [26] have given the expressions for the electromagnetic field of a 
fast charge passing through a slab. We shall briefly repeat the method of obtaining of these 
expressions. 

Since the problem is homogeneous both in time and directions in each domain z < 0, 
< z < a, and z > a normal to the charge velocity, it is convenient to represent all field 
components as Fourier integrals over time and transversal coordinates r = (x,y). Then the 
Fourier component of electric field is obtained from the Maxwell equations: 
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where n = u/u, k = (k x ,k y ), po(k,uj,z) is the Fourier component of the charge density of the 
particle 

Po(k '^' z)= (2^ eXP (v)- (2) 
The Fourier component of the magnetic field is expressed through E(k, lo, z) as follows: 

B(k, u,z) = - {-i [V x E(k, u, z)\ + [k x E(k, u, z)]} . (3) 
ijj 

The total solution of (1) for a charge density (2) is a sum of solutions to homogeneous and 
inhomogeneous equations. While the first equation describes the radiation field, the second 
equation describes the particle field proper in a medium with local dielectric function e(tu). 
Also, equation (1) must be solved for each domain inside and outside the slab, and therefore 
the solutions are joined using the boundary conditions (equality of normal induction components 
and transverse electric field components on the boundary) 



E z (k,u,-0) = e(tj)E z (k,Lo,+0), e(u))E z (k,w,a-0) = E z (k,u,a + 0), (4) 
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k-E(k,w,-0) = k-E(k,w,+0), k-E(k,w,a-0) = k • E(k,w, o + 0). 
Taking account of these conditions, the following system of relations is obtained: 



E(k,w,z) = E (1) (k, uj) exp [i—z] + 77— ^fj a[ ) (kn + k^-ri) exp (-i-nz) , (5) 
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D(k, uj) = (r 2 + et\) 2 exp ^-i^r 2 a^ - (r 2 - en) 2 exp ^i^r 2 a^j (16) 

and = u/c, r )~ 2 = 1 — (3 2 . The functions a 2 ^ and a 2 + ^ are expressed through a[ ^ and 
and are not explicitly given here. They may be obtained from the matching conditions for the 
normal component of the electric induction on the surfaces z = and z = a. 

The first terms in (6) and (8) describe the Coulomb field of the particle. The first term in 
(7) describes the particle field in an unbounded medium characterized by the dielectric function 
e(uj). The field is identical with a Cherenkov radiation electric field in the frequency range 
(3 2 e{lo) > 1. All other terms are due to existence of boundaries. Particularly, they describe 
the transition radiation in the backward (second term in (6)) and forward (second term in (8)) 
directions [26]. 

3 The wake charge evaluation 

In this section we shall consider the volume charge induced by a moving particle in a slab (the 
so-called wake charge). To evaluate the wake charge, Poisson equation is used 

p v = (1/4tt)VE - po (17) 

in which the pQ = qS(r)S(£) is the charge density of a test particle, £ = z — ut, E is the electric 
field in the slab which is determined by the inverse Fourier transformation of (7). Since the 
divergence of the second term in (7) (of the radiation field) is zero, the wake-charge density is 
determined only by the first term in (7). 

Using this term in the relation (18) we obtain 
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ft = ^-iW/ +0O <ia,expto)l^M, (18) 

2lXU J-oo \ u J e[u) 

where 8{x) is a Dirac function. Since the dielectric function of the medium has poles only in 
the lower u half-plane [6], no induced charge exists in front (£ > 0) of the particle. Note that 
the relation (19) may be also obtained from the expression for electrostatic potential created 
by the particle in unbounded medium described by a dielectric function e(u) [6]. 

Evaluating an integral of expression (19) over the volume we obtain the wake charge, induced 
by the particle, moving in the slab: 



Q v (t) = -qMt)-Hr)], (19) 



where 



*(*) = -W + °°^exp (--,*) i^M= (20) 

= "K 1 _ + m { 1 _ h ~ ? exp [Aj cos ^ + Sj sin {ujjt)] } 

and r = t — a/u, eq is the static dielectric constant of the medium, 9(t) is the Heaviside unit step 
function (with 6(0) = ^), the symbol P denotes the principal value of the integral, ±ujj — ivj 
are the solutions of the equation e(u) = (vj > 0), while coefficients Aj and Bj are 

A j = 2Re\- — \- — 1, B j =21m\- — i- — 1. (21) 

I [ujj — iv j) e (Uj — iv j) \ I (Uj — iv j) e (Uj — iVj) \ 

Here the prime denotes differentiation with respect to the argument. The summation in (21) is 

carried over all zeros of the dielectric function. 

Analytic properties of the dielectric function [6] and the residue theorem were used in 

evaluation of 3>(i). 

4 Calculation of the induced surface charges 

The surface- induced charge density is related to the discontinuity of the electric field z-component. 
Expressions (6)- (8) give 

//■ + 0O 
oPk / duai(\t, w)exp[i(kr — ut)], (22) 
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and indices i = 0, a refer to the first and second boundaries. 

The total induced charge is obtained by integration of expressions (23)- (25): 

Qis{t) = / d 2 rai(r,t) = (2tt) 2 f doj<Ji{uj) exp(— iujt), 

J J — oo 



(24) 



(25) 



where (Ji(uj) = <Tj(k = 0,uj). 

Let us consider the dynamics of induced charges for the dielectric slab, the dielectric function 
of which, as it is known, has no singularity in the static limit, when to — > [6]. It is clear 
from expressions (12) and (13) that the functions a[ -* and in expressions (24) and (25) are 
proportional to k 2 , and calculating the function <Ji(uj) for the dielectric slab the above-mentioned 
functions tend to zero. Thus, the charges being induced on the surfaces of the dielectric slab 
are determined only by the first terms in the expressions (24) and (25), that is by the electric 
fields which are created by the particle in an unbounded dielectric with the dielectric function 
e{lo) and in the vacuum when the dielectric is absent. 

Using the known Sokhotsky-Plemel relations [6] for quantities Oi(uj) in (26) we have 
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The following relations are obtained from expressions (21), (26)-(28) 



(26) 
(27) 



Q os (t) = q$(t), Q as {t) = -<?$(r). (28) 

Thus, expressions (21) and (29) make it possible to obtain the total surface charge as soon as 
the zeros of the function e(uj) are known. One may easily verify that, at any time, the net 
induced charge (a sum of Q os , Q as and Q v ) in the slab is zero. 

The following interpretation for expressions (20), (21) and (29) may be given. When the 
particle approaches the slab surface from a vacuum (t < 0) we have Q os = —(q/2)(l — 1/eq), 
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Qas = ((//2)(1 — 1/eo)) and Qv = 0. Note that the first boundary of the slab is charged 
oppositely to the second boundary. 

While the particle moves inside the medium (0 < t < a/u), the first boundary charge 
oscillates and decreases: 

Qos(t) = I (l - ^) - q Yl ex P i A J cos + B i sin • ( 29 ) 

Meanwhile the volume charge is increased: 

The second boundary charge remains unchanged. 
For t > a/u we have 

Qas(*) = -| ^1 - j^j + g^exp(-iyr) [A,- cos (w/r) + sin (ujjt)] , (31) 



Qu(i) = 1^2 ex P (~ v jt) \Aj cos + sin (ujjt)] — (32) 
i 

~~ 1 ex P ( — ^i 1 ") IA? cos i u} j T ) + s i n ( 0J j T )] • 
i 

The charge on the first boundary is given in this case by expression (30). 

After the particle crosses the second boundary (t > a/u), the charge on the first boundary 
decreases to its lower limit (q/2)(l — 1/sq). On the second boundary the charge value increases 
and attains its maximum —(q/2)(l — 1/eq) when t 3> a/u. The wake charge in the volume 
becomes equal to zero. 

Thus it follows that after the particle crosses the first boundary, the surface charge is trans- 
formed into the wake charge. The latter is transformed again into the surface charge after the 
particle crosses the second boundary. 

We apply the results represented by expressions (20), (21) and (29) for the induced charges 
in the model of diatomic cubic ionic crystal or polar semiconductor, whose dielectric function 
is given by [27] 



8 



e(lo) = e, 



lo 2 t — lo 2 — ivuj 



(33) 



oo 



lo t — uj 2 — ivlo 



In this expression is the optical frequency dielectric function, lot and lot are the frequencies of 
the longitudinal and transverse-optical vibration modes of infinite wavelength, v is the damping 
rate, which we assume to be small (v <C lol), and £q is the static dielectric function, which enters 
the theory through the Lyddane- Sachs- Teller relation, lo t Eq = lo^e^. Expression (34) implies 
that its root e(lo) = has the following form: 



where O, 2 = lo\ — v 2 /4. Substituting (34) and (35) into expression (22), we find the coefficients 
determining the net induced charges: 



Figure 1 depicts the time dependence of Q os (t), Q as (t), and Q v (t) for the slab of LiF. The 
following values of the parameters were taken for numerical calculations: = 1.96 and Eq = 
9.01 [27], v/lol = 0.2, and Qa/u = 15. We see that as the particle crosses the boundary, the 
surface and wake charges oscillate with a frequency Q, although the net induced charge remains 
equal to zero. 

5 Conclusion 

Let us briefly discuss the conditions in which the processes taking place at the boundaries of the 
slab can be considered independent and the boundary can be interpreted as that of a half-space, 
as was done by Gorbunov et al [19]. 

We see from expressions (20), (21) and (29) that if the condition a < u/vj is met, the charge 
— q(l — 1/eo) has no time to transform into the wake charge before the particle reaches the second 
boundary. For this reason the transformation of a surface charge into a wake charge and the 
transformation of the latter into a surface charge at the second boundary are interrelated. When 
the particle crosses the second boundary of the slab, near the boundary it excites electric-field 
oscillations [3, 20, 26] whose phase is related to that of the oscillations of the electric field near 



LOj — IV j = — 




(34) 




(35) 
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the first boundary. For the fields as the boundaries to be completely independent, a must exceed 
u/vj. In this case not only the amplitudes but also the phases of oscillations of the electric 
fields at the boundaries are independent. 

We would like to end our consideration with the following note. The densities of induced 
surface charges (expressions (23)- (25)) and its corresponding surface current densities may be 
considered as the sources of transition radiation [27]. As follows from (24) and (25), the densi- 
ties of induced surface charges are determined by both electrostatic and electromagnetic fields 
created by the moving particle. The total induced surface charges, determined by (21) and 
(29), do not depend on the velocity of light and, therefore, are determined only by electrostatic 
fields. Thus, electromagnetic fields do not contribute into the total induced surface charge. It 
seems that relativistic effects apparently reveal themselves in the case when the thickness of the 
medium traversed is limited. 
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Figure Caption 



Fig.l. The dynamics of the induced charge at the front boundary (full curve), the rear 
boundary (broken curve), and in the volume (dotted curve) for the slab of LiF. The following 
values of the parameters were taken for numerical calculations: = 1.96, £0 = 9.01, v/lol = 
0.2, and Qa/u = 15. 
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